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ABSTRACT 

We study merger histories of dark-matter haloes in a suite of TV-body simulations that span dif- 
ferent cosmological models. The simulated cases include the up-to-date WMAP5 cosmology 
and other test cases based on the Einstein-deSitter cosmology with different power spectra. 
We provide a robust fitting function for the conditional mass function (CMF) of progenitor 
haloes of a given halo. This fit is valid for the different cosmological models and for different 
halo masses and redshifts, and it is a significant improvement over earlier estimates. Based 
on this fit, we develop a simple and accurate technique for transforming the merger history of 
a given simulated halo into haloes of different mass, redshift and cosmology. Other statistics 
such as main-progenitor history and merger rates are accurately transformed as well. This 
method can serve as a useful tool for studying galaxy formation. It is less sensitive to the 
low accuracy of the fit at small time-steps, and it can thus replace the more elaborate task 
of construction Monte-Carlo realizations. As an alternative approach, we confirm the earlier 
finding by Neistein & Dekel that the main-progenitor follows a log-normal distribution. This 
property of merger trees allows us to better capture their behaviour as a function of time and 
descendant mass, but a broader suite of simulations is required for evaluating the dependence 
of the log-normal parameters on the cosmological model. 

Key words: cosmology: theory — dark matter — galaxies: haloes — galaxies: formation — 
gravitation 



1 INTRODUCTION 

The growth of dark-matter haloes through merging and accretion 
is driving the formation and evolution of galaxies. An accurate 
theoretical prediction for the way haloes grow is thus a key ele- 
ment in the effort to develop a theoretical understanding of the pro- 
cesses associated with galaxy formation, including star formation, 
the growth of black holes in galaxy centers and their appearance as 
quasars. 

The conditional mass function (hereafter CMF) has been an 
important tool in quantifying the growth of haloes. It is defined at 
a given time as the average number of progenitors that will merge 
into a descendant halo at a later time. The CMF was introduced 
in the context of the Extended Press-Schechter formalism based on 



excursion- set theory (hereafter EPS. lBond et al.ll99lMBower|1991l ; 
lLacev&Coiej|l993l) . Recent theoretical predictions use a variant 
of the EPS formalism in which the sphe rical collapse model is 
replaced by an ellipsoidal collapse m odel JSheth & Tormenll2002l : 
iMoreno et al.l2008l : IZhang et aU2008h . An em pirical fit to the CM F 
based on A^-body simulation was presented bv lCole et al.l J2008l> . 

Although the excursion-set models that employ ellipsoidal 
collapse are successful in reproducing the overall (unconditional) 
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halo mass function as derived from A^-body simulations these 
models do not provide accurate CMF predictions. Most of the cur- 
rent analytic predictions of the CMF deviate from the results of N- 
body simulations by a multiplicative factor of a few, es pecially for 
the nu mber of massive progenitors. The empirical fit of ICole et al.l 
d2008h does better, but its inacc uracies may reach the level of 50%. 
As shown bv lCole et al] J2008h . these inaccuracies are valid when- 
ever the standard variables are used, and they cannot be reduced 
by a different functional fit. In ad dition, this fit was calibrated 
against the Millennium simulation JSpringel et alll2005h that was 
based on the cosmological parameters from the first-year data of 
the Wilkinson Microwave Anisotropy Probe (WMAP), it should be 
interesting to evaluate the CMF for the more up-to-date cosmolog- 
ical parameters that emerge from the fifth-year WMAP data (e.g. 
lKomatsuetaI]|2009l) . 

The main goal of this paper is to provide a more accurate and 
more robust empirical description of the CMF as measured from 
Af-body simulations. We work out the possible scaling laws which 
can be applied to the CMF in order to capture its detailed proper- 
ties over a large range of cosmological models, halo masses and 
redshifts. The result is a fitting function that offers a significant im- 
provement in accuracy over previous studies. The empirical fit can 
help us distinguish between different analytical models of structure 
formation, and can guide us to new improved versions of them. 
Such a fit can be used for generating Monte-Carlo merger trees, 
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which will accurately reproduce the results of TV-body simulations. 
The effect of cosmology, environment density and different dark 
energy model s can be studied, re lating these factors to haloes and 
galaxies (e.g. iMaccio et al.l|2008i ). These issues and additional ap- 
plications are discussed below in more detail. 

Currently, there is still some freedom in the definitions of halo 
mass and merger trees. Accurate CMF fit may help us to test dif- 
ferent algorithms for constructing merger trees, including assign- 
ing masses to haloes and relating progenitors to their descendants. 
It will allow us to quantify the effects of different construction 
schemes for merger trees, and especially to identify the scheme that 
results in the simplest and most universal CMF. 

A significant progress has been made in quantifying the over- 
all halo mass function, namely the abundance of haloes of a given 
mass a t a given redshift. The analytical model of Sh eth & Tormeri 
i 19991) offers a significant improvement over the classic estimate of 
IPress & Schechteii J 19741) . Furthermore, the growing volume and 
dynamical range of iV-body simulations allow an accurate mea- 
surement of the halo ma ss function with small sampling scatter 
and cosmic variance (e.g. Ijenkins et al]|200ll : IWarren et al]|2006l ; 
iReed et alj2007h . Still, there are significant deviations between the 
analytic approximations and the simulations, which can increase as 
a function of redsh ift and reach a few tens of percents at z = 2.5 
dTinker et alj|200^) . For one thing, an accurate knowledge of the 
CMF can provide a better estimate of the halo mass function and 
its evolution at high redshift. However, the more interesting fea- 
ture of the CMF is that it involves much more information than the 
halo mass function, and it can still be derived with a manageable 
statistical uncertainty. 

Using the CMF, we present a new method for scaling a given 
set of merger trees into a different cosmological model, mass res- 
olution or redshift. Although this technique requires an existing 
database of merger trees extracted from a given A-body simu- 
lation, it offers an easy transformation to the desired cosmology, 
mass and redshift. This method has certain advantages over gener- 
ating Monte-Carlo trees based on excursion-set models. For exam- 
ple, it bypasses the inaccuracies associated with the non-Markov 
features of the trees, associated wi t h correlations between ad jacent 
time-steps JSheth & Tormell2004lN~eistein & De kel 2008a). This 
method can be extended to include halo substructure and halo loca- 
tions. Such a method can be very useful for semi-analytical mod- 
eling of galaxy formation, black hole growth, and dwarf galaxy as- 
sembly. 

The common line of thought in quantifying merger trees starts 
with the CMF. This is the basic prediction of the EPS formalism, 
and it includes substantial information on the statistical properties 
of the tree. The known CMF allows derivations of other statis- 
tics such as the main-progenit or history and merger rates. How- 
ever. iNeistein & Dekel d2008al . hereafter ND08a) have shown that 
an alternative path of reasoning may also be useful, where the 
main-progenitor history is derived first and a full construction of 
merger trees follows. They found in A-body simulation that the 
main-progenitor statistics is very regular and follows a log-normal 
distribution of an appropriate mass variable, the variance S of the 
smoothed density field. One of our purposes in the current paper 
is to verify this log-normal behaviour in a range of cosmological 
models. We hope that this will bring more insight into the physical 
origin of the main-progenitor statistics. 

It should be noted that the CMF does not include all the in- 
formation needed for describing merger trees. There are in general 
many different subsets of trees that can accurately fit a given CMF. 
Nonetheless, we will show in this work that the scaling laws of 



Table 1. A summary of the ./V-body simulations used in this work, all with 
flat cosmology and with a Hubble constant of 72 km/s/Mpc. Particle mass 
is in units of h —1 MQ, box size is in Mpc. Each simulation follows the 
evolution 350' 5 particles. 



Name 






Particle 


mass 


Box 


wmap5 


0.258 


0.796 


4.54 X 


10 s 


90 


lcdm 


0.258 


0.915 


4.54 X 


10 8 


90 


scdml 


1.0 


0.77 


1.76 x 


10 9 


90 


scdm2 (n = — 2) 


1.0 


0.8 


1.76 x 


10 9 


90 



the CMF can provide a good estimate for the full statistics of the 
merger trees, including main-progenitor histories and merger rates. 

This paper is organized as follows. In ij2]we describe the suite 
of /A-body simulations used and the way merger trees are con- 
structed. Section [3] is devoted to the conditional mass function, 
where we study its scaling properties and provide an accurate fitting 
function. In Sj4]we discuss an alternative description, using the log- 
normal nature of the main-progenitor history. A simple prescrip- 
tion on how to scale a given simulation is developed in S|5] We 
summarize the results and discuss them in $6] Additional statistical 
properties of merger trees which are not critical for the body of the 
paper are added in Appendix [A] Throughout the paper we use log 
for log 10 and In for natural logarithm. 



2 THE SIMULATIONS 

All simulations have been performed with PKD GRAV, a tree code 
written by Joachim Stadel and Thomas Quinn dStadelll200ll) . The 
code uses spline kernel softening, for which the forces become 
completely Newtonian at 2 softening lengths. Individual time steps 
for each particle are chosen proportional to the square root of the 
softening length, e, over the acceleration, a: Ati = if^/ejai. 
Throughout, we set r\ = 0.2, and we keep the value of the soft- 
ening length constant in co-moving coordinates during each run 
(e=1.62 h, _1 kpc). Forces are computed using terms up to hexade- 
capole order and a node-opening angle 9 which we change from 
0.55 initially to 0.7 at z = 2. This allows a higher force accu- 
racy when the mass distribution is nearly smooth and the relative 
force errors can be larg e. The initial condit ions are generated with 
the GRAFIC2 package (Bertschinger 2001). The starting redshifts 
Zi are set to the time when the standard deviation of the smallest 
density fluctuations resolved within the simulation box reaches 0.2 
(the smallest scale resolved within the initial conditions is defined 
as twice the intra-particle distance). For each simulation we stored 
more than 100 outputs from redshift 10 to redshift zero, in order to 
construct detailed merger trees. The parameters of the simulations 
used in this work are describe in TableQ] 

In all sim ulations, dark matt er haloes are identified using the 
FOF algorithm jDavis et"al ] [l98i with linking length of 0.2 times 
the mean interparticle separation. Only haloes which include more 
than 20 particles are saved for further processing. For constructing 
merger trees, we started marking all the particles within the virial 
radius of a given haloes at z = and we tracked them back to the 
previous output time. We then make a list of all haloes at that earlier 
output time containing marked particles, recording the number of 
marked particles contained in each one. In addition we record the 
number of particles that are not in any halo in the previous output 
time and we consider them as smoothly acc reted. 

We used the two criteria suggested in IWechsleretai] J2002t) 
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Figure 1. The power spectra of the density fluctuation fields in the simula- 
tions. Shown is S(M) = cr 2 (M), the variance of the initial mean density 
in spheres that encomp ass mass M, linearl y extrapolated to z = 0. The 
Millennium simulation (Springel et alj2005l) is shown for reference. 



for halo 1 at one output time to be labeled a "progenitor" of halo 
2 at the subsequent output time. In our language halo 2 will then 
be labeled as a "descendant" of halo 1 if i) more than 50% of the 
particles in halo 1 end up in halo 2 or if ii) more than 75% of halo 

1 particles that end up in any halo at time step 2, do end up in halo 

2 (this second criterion is mainly relevant during major mergers). 
Thus a halo can have only one descendant but there is no limit to 
the number of progenitors. 

We found evidence for the so-called 'backsplash' subhalo 
population (e.g. iKnebe et al, I l2008l) . These haloes have orbits that 
brought them inside the virial radius of their host at some earlier 
time, but without them been really accreted (i.e. they managed to 
come out from their host dark matter halo). We decided to treat 
them in two different ways according to their final fate: i) If af- 
ter have been inside the main halo, the backsplash halo survives as 
isolated halo till the present time, than it is removed from the pro- 
genitor list of the parent halo (i.e. it is removed from the merger 
tree); ii) if the halo is accreted again (in a definite way) by the main 
halo at a later time step then it is considered as accreted the first 
time it entered the main halo. We found that back splash haloes are 
roughly 8% of the total progenitor number but they only marginally 
contribute (less than 2%) to the final halo mass. 

According to the Extended Press-Schechter formalism the sta- 
tistical properties of merger trees are fixed by the density field 
at early times when perturbations grow linearly. All the statisti- 
cal properties of this field are described by its variance, S(M) = 
a 2 (M). Specifically, S(M) is the variance of the density field, 
smoothed with a spherical top hat filter i n space, and linearly ex- 
trapolated to z = (for more details, see lLacev~ & Cole 1993b- m 
fig [Il we plot S(M) for all the simulations used in this work. This 
comparison shows the predicted similarity between merger trees of 
different simulations: we have two simulations with very similar 
S(M) but different £7 m (lcdm & scdm2). One simulation has a low 
value of erg (wmap5), and the scdml simulation has a very different 
shape of S(M). 

In order to avoid deviations of S(M) due to the small box size 
and cosmic variance we measured it directly from the initial condi- 
tion density field used for each simulation. The values of as given 



in table [TJ are those obtained by this calibration method. For the 
scale free simulation (scdm2) the box size is important as it limits 
contributions from large scales to the field variance (see a detailed 
discussion in lSmith et alj|2003l) . This effect bends S(M) slightly 
at the high mass end, in agreement with the theoretical prediction. 



3 CONDITIONAL MASS FUNCTION 

This section is devoted to a comparison of the conditional mass 
function (CMF) as extracted from our suite of iV-body simulations. 
The scaling laws of the CMF are important as they highly constrain 
the behaviour of the full statistics of merger trees. In Appendix lAl 
we show other statistical properties of trees, and demonstrate that 
they follow approximately the properties of the CMF found here. 

According to the Extended Press-Schechter for malism (EPS, 
iBond et aljfl99lt lBowerlll99ll : lLacev & Colej[T993h . the average 
number of progenitors in the mass interval [M, M + dM], which 
will merge into a descendant halo Mo after a timestep Aui, is given 
by 

^(M|M 0) A W )dM = ^/(AS,Au,) £L dM . (1) 



Here ui = 8 c (z)/D(z), where S c (z) — 1.68 with a weak de- 
pendence on z, and D(z) is the cosmological linear growth rate. 
Au; = uj — ujo where the progenitors are identified at ui(z) and 
the descendant halo at loq. We refer the reader to the appendix of 
ND08a for a detailed summary on how to compute ui, and a simple 
fitting function. The quantity / dS describes the fraction of mass 
of the descendant halo that is in progenitors of mass in the range 
[S, S + dS]. According to EPS 



/ ps (AS,Au;) = 



Au 



^ (AS)'/ 2 



exp 



2AS 



(2) 



where AS = S(M) — S(Mo), and / ps is the specific solution of / 
given by the EPS formalism. In what follows we will use the name 
'CMF' to designate dN/dM only. 

In the language of the EPS formalism / ps is 'universal' - it 
does not explicitly depend on the descendant halo mass Mo or 
on the background cosmology. These parameters appear only in 
the transformation of AS and Au; to mass and time (or redshift). 
The limited success of the EPS formalism in reproducing the iV- 
body simulation merger trees indicates that the function / ps is a 
limited approximation, bu t it is possible in princi ple t hat / PB has 



univer sal properties. Both ISheth & Tormenl ( 120021) and lCole et al .1 
d2008l) found that / as derived from A r -body simulations shows 
non-negligible deviations from universality. However, these stud- 
ies did not explore in detail the breakdown of universality and did 
not test other models f or / against iV-body simulations. For exam- 
ple, ICole et al I d2008l) have studied / only for descendant haloes 
identified at z — and showed deviations of ~ 50% for various 
values of Auj and halo masses. 

The universality of / ps is based on the two main theoretical el- 
ements of the EPS formalism, namely the dynamics of the spherical 
collapse model and the statistical properties of the initial density 
fluctuation field when smoothed with a top-hat filter in fc-space. 
Variants of these elements have been attempted, incl uding ellip- 
soida l collapse or smoothing filters of different shapes (|Bond et all 
1 199 it ISheth & Tormenll2002l : IZentnerll2007l : lDesiacquesll2008l) . In 
these models / may in principle depend on the descendant halo 
mass and thus become non-univers al, but this is not necessarily so. 
For example, iMoreno et al. I J2008h showed that when ellipsoidal 
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collapse is assumed, / remains universal if the variables S(M) 
and Auj are normalized by S(Mo) and y S(Mo) respectively. In 
its most general form, / may depend on So = S (Mo), wo and the 
specific cosmology used. We will show below that neglecting the 
dependence on wo and on cosmology may not be too harmful, but 
the dependence on So must be properly taken into account for an 
accurate description of the CMF. 



3.1 Self-similarity in time 

According to eq.[T] the CMF as predicted by EPS depends only on 
Auj and not on the redshift zo where the descendant halo is defined. 
This self-similarity implies that merger trees extracted at different 
redshifts are self- similar when using lo as the time variable. In- 
deed, ND08a and lGeneletail (2008) have verified this behaviour 
usin g merger trees extrac ted from the large Millennium simula- 
tion JSpringel et al.l 120050 . It was shown to work well for main- 
progenitor histories a nd merger rates, with s catter of a few percent 
up to z ~ 5 (see also lFakhouri & mH|2008| . who found deviations 
from this symmetry using a different definition of the halo mass). 

In fig. [2] we show the level of self-similarity in time in our 
four iV-body simulations. The scdml and scdm2 simulations show 
small variations in the CMF between different values of zo, con- 
sistent with the sampling noise. The wmap5 and lcdm cosmologies 
show larger deviations, up to a factor of two for small progenitors 
and for the time-step Auj — 0.5. However, the data do not show 
a monotonic trend with zo, hinting that our small-number statistics 
may contribute to the scatter. For example, the number of descen- 
dant haloes within a mass range 10 12 < M < 10 13 /l _1 M 
in the wmap5 cosmology is roughly (500,400,100,100) for zo = 
(0, 1, 2, 3) respectively. It is encouraging that better results were 
obtained for the Millennium simulation (as described above) and 
for large time-steps. Nonetheless, part of the deviation from self- 
symmetry may be related to the way haloes are defined using FOF 
linking length and the variation with redshift implied by this defi- 
nition. As explained in section[2] we try to correct the merger trees 
for 'backsplash' haloes. This correction may introduce some addi- 
tional asymmetry between zo = and higher redshifts. 

The fact that merger trees remain similar when the descen- 
dant time is varied from z = 0toz = 3isa non-trivial result. 
The different values of the cosmological parameters f2 m and Oa 
at high redshift imply that the merger trees for a higher zo evolve 
in a different cosmological background. Thus, the dependence of 
merger trees on the cosmological parameters Sl m and Qa should 
be all folded into the time variable ui = S c /D. In fig. [3] we test 
this hypothesis explicitly by comparing the CMF from the scdm2 
and lcdm simulations. These simulations have a similar shape of 
S(M) (see fig. [TJ but very different fi m values. The agreement is 
within the Poisson sampling error-bars, demonstrating that uj can 
scale properly different expansion histories of the Universe. 



3.2 Different power-spectrum 

Following the EPS formalism (eq. |2j the fraction of mass in- 
side progenitors, f ps (AS, Auj), depends only on Auj and AS. 
As mentioned above, one could assume less universal form, 
/(AS, Alu\So) that will enable accurate description of the CMF 
for various cosmological models. In order to test this we plot / in 
fig. [4] for two sets of power-spectra. The results of the first com- 
parison (left panels) between lcdm and wmap5 simulations is very 
good, showing no significant deviation in / between these mod- 
els. This is obtained when the same values of Mo are selected in 



scdm2 




lcdm 



Figure 2. Self-similarity in time of the progenitor conditional mass func- 
tion (CMF). Each panel shows the CMF as measured from a different IV- 
body simulation as indicated. The descendant halo is identified at differ- 
ent times corresponding to redshifts zo = 0, 1, 2, 3 (solid, dashed, dotted- 
dashed, and dotted lines respectively). The progenitors are identified at two 
lookback times corresponding to Auj = 0.5,3 (blue and red lines re- 
spectively). The distribution of descendant masses is selected in the range 
10 12 ^ Mo $C 10 13 Mq such that it produces the same distribution 
of Mo in the different cases of zo ■ In each case, the four curves show only 
a weak dependence on the descendant time zo ■ 
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Figure 3. Self-similarity of the CMF with respect to the cosmological ex- 
pansion history. The two sets of CMF are computed for zo = 0, for 
descendant mass in the range 10 12 $C Mo $J 10 13 h~ 1 MQ, and for 
Au; = 0.3, 3 (blue and red lines respectively). The solid and dashed lines 
refer to the scdm2 and lcdm simulations respectively. The distribution of 
descendant masses at zo is matched in order to produce the same distribu- 
tion of S(Mo) values, and the average descendant halo mass is indicated in 
units of /i _1 Mq. Note that we plot the number of progenitors per unit of 
S, which is given by eq.[T]times \dM/dS\. The CMF shows only a weak 
dependence on cosmology. 




Figure 4. Fraction of mass in progenitors for different power spectra and 
cosmologies. The lines of two types in each panel compare results from 
two different simulations as indicated. The time steps are Aw = 0.3, 2 
and the descendant haloes are selected at zo = 0. In the left and right 
panels the descendant halo masses in wmap5 were selected in the range 
10 13 < M < 10 14 h- 1 M e and 10 12 < M < 10 13 /i" 1 M respec- 
tively. In lcdm and scdml the descendant haloes were chosen to match the 
distribution of Mq or So in wmap5, as indicated in each panel. 



each cosmology, or the same So = S(Mo). The difference be- 
tween these two selection criteria is negligible, so we actually sam- 
ple f(AS, Aw | So) for a small range in So- 

The second comparison (right panels) shown in fig. [4] is for 
wmap5 and scdml cosmologies. As seen in fig.fT]both S(M) and 
its derivative dS/d-M vary significantly between these two simu- 
lations. This is being translated into a large discrepancy in / when 
the same values of Mo are chosen from both simulations. How- 
ever, selecting a sample with the same So values leaves / invariant, 
proving that / = /(AS, Aw | So) (at least for our limited set of 
data). For a given So the values of Mo in both simulations differ 
by a factor of ~10, limiting the dynamical range for which we can 
check this / scaling. 

The dependence of / on So deserves a more careful test, prob- 
ably with a larger set of TV-body simulations, spanning a larger 
range of power-spectra. It is important to verify that such a depen- 
dence is not related to the shape of S(M), as was found here. If 
this is true, it implies that a modification to the normalization of the 
power-spectrum (as) which changes So, will induce a non-trivial 
change to the CMF. Such a dependence on as is different from the 
simple scaling of the shape Aw/VAS that is used by EPS (see 
below). 



3.3 Global fitting function 

The EPS formalism suggests that the mass (AS) and time (Aw) 
variables can be scaled into a new variable v = Aw/V AS. In 
terms of v, the function f(y) v dS/df should be universal. As 
shown above some dependence on So is required, so we can try 
to use a function of the shape f(v\So). This is tested in fig. [5] 
where we plot f(u\So) for a fixed value of So and for various 
Aw. The clear deviations from a unique line show that such a for- 
mulation is too simple and incompatible with the results of N- 
body simulations. This figure is very similar to the one given in 



Figure 5. The limitation of using a fitting function of the form f(v\So). 
The results form the scdm2 simulation are plotted as histograms. They are 
extracted at z,o = with three different time steps Aw = 0.5, 1, 3 and for 
an average halo mass of 3 X 10 12 h~ 1 Mq. The global fit of eq.[5]for the 
three time steps is shown as smoothed sol id curves. The EP S prediction is 
plotted as a thin dashed line and the fit of ICole et alj J2008h is shown as a 
dotted-dashed line. Like any other function of the type f(v), they both fail 
to reproduce the simulation results. The lognormal global fit to Pi (eqs.[5] 
&|6) is plotted as small circles (only for M > Afo/2, where it should be 
identical to the CMF, see section[4). 



ISheth & Tormenl d2002l) and ICole et al.1 ( 120081) . and agrees with 
these previous results. Thus, we are forced to use a more general 
form / = /(AS, Aw|S ). 

We found that the fitting function below can fit the CMF 
from our set of TV-body simulations for all time-steps larger than 
Aw ~ 0.5 and all So values. We have looked for the simplest pos- 
sible function that can fit the data, which is still similar to the EPS 
original function. The function we adopt is: 



/(AS,Aw|S ) = 
Aw 2 

exp 



'AS 3 



Aw 



-d 



exp 



Aw 2 



AS 



+ 



(3) 



AS 24 



where the parameters explicitly depend on So in the following way 

0.215 + 0.0037 So 
0.414 + 0.0013 So 
0.0746 + 0.0382 S 
0.257 + 0.0568 S 



(4) 



= -0.0141 So + 0.0056 So 



The added term is not a natural feature of the standard EPS for- 
malism, and it emphasizes the need for additional physics, beyond 
the Sheth-Tormen generalization from spherical to ellipsoidal col- 
lapse model. The failure of functions that depend solely on v in re- 
producing the detailed evolution of TV-body simulations may have 
interesting theoretical implications. We also note that our fitting 
function is not necessarily unique — there may exist other fitting 
functions of comparable quality. However, it is likely to be suffi- 
cient for most applications related to galaxy formation, based on 
the tests performed here over wide ranges of halo masses, redshifts 
and cosmological models. 

In fig. [6] we compare our fitting function to the simulation 



6 E. Neistein, A. V. Maccid & A. Dekel 



0.1 



0.01 



0.1 



0.01 



I 1 1 1 

k S =5.2 ■ 
J \ Au=0.5 


i i i i 

W So=8.9 - 
i% Aw=0.5 

4 1 '"I""- ^I^t^r^ 


fit 

-ft EPS So=1.5 ■ 

0. ST A.' 0.5 

J \ - - Cole 

. ^ S ' Ti; ^ ; 0sLTl= s 


A S =2.3 ■ 
\\ Aw=0.5 



0246 0246 

AS 




Figure 6. A global fit to /(AS, Au>\So) as defined in eqs. [3]and|4] Each 
panel refers to a given set of So and Auj values as indicated. The simulation 
data (histograms) were extracted and averaged from all the different simu- 
lations used here at = and at zo = !■ The descendant masses were 
selec ted in log bins of size half dex. Our global fit, the EPS p rediction, the 
fit of lCole et alj J2008h . and the fit of lSheth & Tormenl j2Q02b are shown as 
solid, dashed, dashed-dotted an d dotted lines res pectively. Note that both 
the EPS prediction and the fit of lCole e t al.| fcOOgl) are independent of So . 



data, showing that it accurately reproduces the trends found. The 
agreement between the fit and data is at the level of the statisti- 
cal noise. The dependence on So is mainly seen at low AS (corre- 
sponding to massive p rogenitors) . We al so show the EPS prediction 



20081). bot h indep endent on So- 
Cole et alj dwS) for intermedi- 



and the fit derived bv lCole et alj 
Evidently, our fit agrees well with 
ate and massive haloes (although deviations for massive progeni- 
tors are highlighted in fig. [5}. This is encouraging because the fits 
are based on different simulations and different merger-tree con- 
struction schemes. Ho wever, our fit bre aks the symmetry of using 
only v as was done bv lCole et al] j2008h . in this way it can capture 
the behaviour of the data for low mass haloes, and across various 
time-steps. The accuracy for different Acu can be specifically seen 
in fig. |H 

The integral of / can be computed analytically only for a = 2 
(So ~ 2.5), for the whole halo mass range we computed the in- 
tegral numerically, yielding 0.75 up to 0.9, depending on the de- 



scendant mass and on the time-step. This predicts that a substantial 
fraction of the mass is not included in any progenitors but rather 
accreted from a 'smooth' component. 

Throughout this work we mainly discuss the EPS formalism 
in its standard version. However, as mentioned in the introduction, 
versions of this formalism which use the ellipsoidal collapse model 
are presumably more accurate, as is indicated by their ability to 
predict accurate halo m ass functions. We exam ine two such stud- 
ies, the pioneer work o flSheth & Tormenl feOO j) and a more recent 
work bv lMoreno et alj d2008l) . Our interpretation of both studies is 
that the CMF as predicted by the spherical model is comparable 
in its accuracy to the one pre dicted by ellipsoidal co llapse. This 
can be seen in figs. 7 & 8 oflSheth & Tormenl J2002b . More ev- 
idently, figs. 3-5 of iMoreno et alj l l2008h show that the spherical 
model gives better results for low mass descendant haloes, but this 
trend changes for massive haloes, where the ellipsoidal collapse 
is more accurate . We a lso show in fig. [6] the CMF proposed by 
ISheth & Tormenl J2002|) FI It is clear that this model predicts CMFs 
which are significantly higher from the simulation results in the 
regime of massive progenitors and large time-steps. Note that plot- 
ting /(AS, Aoj|So) shows more cle arly this effect than t he usual 
plot as a function of mass ratio (e.g. IMoreno et al]|2008h . In Ap- 
pendix lA3l we provide more results in the form of fig. [6] 



4 MAIN-PROGENITOR HISTORIES 

The "main-progenitor" history of a given merger-tree is con- 
structed by following backward in time the most massive pro- 
genitor in each merger event. This is a useful definition as it al- 
lows us to follow a well defined branch of the tree. In addition, a 
quantitative description of the main-progenitor history highly con- 
strains the full statistics of trees. Properties of main-progenitor 
histories were studie d extensively, both analytically and using 
iV-body simulations jLacev & Cole] 1 19931: iNusse r & Sheth' l 19991; 
Firmani & Avila-R eese 2000; Wechsler et al. 2002; va n den Boschl 
20021; INeistein et alj|2006l ; iLi et alj|2007l ; INeistein & Dekelll2008 j ; 



Zhao et all 20081) . Let us also mention that the main-progenitor is 
not always the most massive progenitor at a given time. 

We define Pi as the probability density to find a main- 
progenitor of mass S(Mi) for a given descendant halo mass Mo, 
and a time-step Alj. As was found by ND08a, Pi can be well fitted 
by a lognormal distribution, 



Pi(ASi|S ,Aw) 



CTpASl V 27T 



exp 



(In ASi — Hp) 



(5) 



Here ASi = S(Mi) — S(Mo), the parameters (a p , n v ) depend on 
So and Au), and as usual So = S(Mo). By definition, the integral 
over Pi(ASi) equals unity for any So and Auj. 

Main-progenitor histories were constructed for all the simula- 
tions used in this work. Following ND08a, we confirm that the log- 
normal distribution fits accurately the simulation results. This was 
tested for all possible ranges of halo masses, redshift, and time- 
steps. As was mentioned by ND08a the fit becomes inaccurate for 
small time steps, typically for Alj < 0.5. Such a behaviour is found 
here as well. One exception to the above, where the lognormal fit 
is somewhat innacurate, is for the scdml simulation at zo — 0. 



We adopt eq. 7 from Shet h & Tormenl |2002) and use n! in the denomi- 
nator of T instead of n. According to J. Moreno (private communications) 
this seems to be a typo. 
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Figure 7. The distribution of main progenitor masses in different cosmo- 
logical models. Results are shown for descendant haloes defined at zo = 1 
and for main progenitors at two different time steps, Alu = 0.5, 3. The de- 
scendant halo mass is in the range 10 12 $5 Mo ?C 4 X 10 12 /i _1 Mq for 
the wmap5 cosmology, and Mo in the other cosmologies were selected to 
match the same distribution of S(Mo). Thin solid lines show the lognormal 
global fit, eq.|6] as computed for the Icdm and scdml cases. 

This may be due to the high sensitivity of AS in this cosmological 
model (dS/dM is typically 7 times larger than other cases). The 
fact that the fit works well for z > 1 data using the same cosmo- 
logical simulation is encouraging. In fig.[7Jwe plot main-progenitor 
distributions for few cases of halo mass, zo and Aw, more exam- 
ples can be found in Appendix [B] The accuracy of the lognormal 
fit as plotted here is typical for the rest of the cases. 

For each cosmology we are able to use a global fit, similar to 
the one suge ested by NDO8£0 The fit approximates the standard 
deviation and average of In ASi by 

cr p = (ai log Sb + a 2 ) log Aw + a 3 log So + «4 (6) 
H P = (6i log So + b 2 ) log Alu + b 3 log So + &4 . 

We provide numerical details for the values of these parameters in 
Appendix [B] For a given So, the evolution with time of o v and n P 
is simply linear in log Aw. In fig. [5] we show that this evolution 
gives a very accurate fit to the relevant part of the CMF (the CMF 
and main-progenitor distribution are identical for M > Mo/2). 
We recall that such time-evolution was obtained for the CMF using 
much more complicated Alu dependence (see eq.[3)- Although both 
approaches discuss only the fitting possibilities, it seems that the 
main-progenitor gives a much simpler way to describe the merger 
history in an accurate way. 

The simple behaviour of Pi for each cosmology as seen in 
eq. [6] and in fig. [5] calls for a more global fit, which predicts the 
parameters of Pi in other cosmological models. We compared Pi 
for our set of simulations and found a high level of self-similarity, 
with the exception of the scdml case. Unfortunately we could not 
find a general law able to combine scdml with the other cosmo- 
logical models. For example, choosing haloes with the same values 
of So from scdml and other simulations does not yield the same 
results. To summarize, our study indicates that the main-progenitor 

2 Here we parameterize the fit according to So in order to stress universal 
features, and not according to Mo as was done in ND08a 



Figure 8. The conditional mass function (CMF), dJV/dM, computed for 
different cosmological models using a transformation of the time-step Alu. 
In each panel we compare results from two different simulations and for the 
same descendant halo mass as indicated (in units of /i -1 Mq). The time 
steps used for the Icdm cosmology are Au) = 0.5, 2. For the wmap5 and 
scdml simulations they are different by a factor of 0.86 and 1.8 respectively. 
We learn that the time transformation allows a reasonable scaling of the 
merger trees, albeit with certain deviations in some cases. 

distribution as extracted from A-body simulations has a univer- 
sal lognormal shape. However, the dependence of the lognormal 
parameters on cosmology is not clear, and there is no theoretical 
explanation to this phenomena ye t. 

Recently, IZhao et al ] 120081) studied the behaviour of main- 
progenitor histories for different cosmological models. They pro- 
vide a fitting procedure for estimating median main-progenitor his- 
tories in any cosmological model. However, we find it hard to esti- 
mate the full distribution Pi (AS) out of their algorithm. 



5 APPLICATIONS 

An accurate fitting function for the CMF includes substantial in- 
formation on halo growth, and it has various interesting applica- 
tions: it can shed more light on the evol ution of the halo ma ss func- 
tion with redshift and cosmology (e.g. iTinker et alj|2008b . it can 
constrain dark energy models, it is useful for generating Monte- 
Carlo merger trees, and for predicting mass-accretion histories of 
the main-progenitor. Here we focus on a new methodology for re- 
scaling a given set of merger trees between different cosmological 
models. We also discuss briefly the possibility of generating Monte- 
Carlo trees. 



5.1 Re-scaling a given merger tree 

There are many existing resources of A-body simulations and re- 
lated merger trees which are publicl y available. The mo st useful 
one being the Millen nium simulation dSpringel et al .120051) with its 
web-based database jLemson et al.l200q) . However, recent changes 
in the observed values of cosmological parameters (especially as) 
make these simulations inaccurate for predicting observable quan- 
tities in our Universe. Here we suggest a new methodology to trans- 
form merger trees into a different cosmological model, halo mass, 
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and redshift. Such transformation can also be useful for enhanc- 
ing the mass resolution of merger trees. There are few benefits for 
this approach over the standard method of generating Monte-Carlo 
trees: (a) it can preserve the non-Markov behaviour of trees (see 
e.g. ND08a) (b) it might be easily extended to handle substructures 
(c) it might be extended to accurate transformation of halo spatial 
locations. 

We define the case r as our reference data from a given 
merger-tree, 



{Mi | Mo, Au>, C}, 



(7) 



The data is defined as a set of progenitors Mj for a given descendant 
halo mass Mo, time-step Au; and a cosmological model C. Our 
target case is defined as a different cosmology C and descendant 
mass Mo 



{Mi\M ,Auj,C} t . 

We are looking for a transformation of the kind 

Mi -> Mi , Aw -> Au 



(8) 



(9) 



that will yield a different set of progenitors, possibly for a different 
time-step, which will be consistent with the target case. Such a con- 
sistency is achieved when the transformed progenitors will yield the 
proper CMF. Note that our global fit was done for /(AS 1 , Au>|So), 
as seen from eq. [Tj the CMF includes two other components, the 
mass ratio Mo/M and the derivative dS/dM. As a result, transfor- 
mations that will yield the appropriate CMF cannot be done using 
the information imprinted in / alone, and the functional depen- 
dence of S on M should be important. 

We start by searching for the most accurate transformation 
in Au; which can compensate for the change in cosmology C, or 
descendant halo mass Mo. This means that we only change Au; 
into Auj while keeping the masses of the descendant and progen- 
itor haloes fixed. In the general case, such a simple transforma- 
tion might not allow us to accurately transform the trees. In prac- 
tice, it might be a good starting point for many relevant cases. Our 
methodology to find time-transformation is as follows: for a given 
reference case we use eq.[3]to compute the reference CMF. We then 
compute the CMF for the target case using various time-steps Auji. 
The value of Au is chosen from the set of Auji in order to give the 
best match between the reference and the target CMF's. 

In fig. [8] we plot the CMF for two sets of cosmologies where 
time transformation is being used. This transformation is relatively 
accurate, although it can introduce a large transformation in time 
(almost a factor of 2 in Au; for one of the cases tested here). The 
original differences between the simulations can be seen in fig. [2] 
for reference. Nonetheless, it seems that the transformation in Au; 
do not provide a uniform accuracy for all descendant masses. In 
fig.[9]we show how changing Au; can compensate for different de- 
scendant halo masses. This indicates that enhancement of mass res- 
olution can be easily obtained. 

Time transformations are limited in accuracy, and can lower 
the time resolution of a merger-tree significantly. In addition, such 
transformations might stretch the non-markov correlations between 
consecutive time-steps in a way that might be different from the 
behaviour of the simulations. Our statistical sample does not allow 
us to explore the last effect in detail. 

A less trivial transformation is needed when the time scaling 
cannot provide an accurate enough solution. In this case we can 
transform the mass of the progenitor haloes in order to yield the 
same CMF. For any progenitor mass M we define M such that the 
integral over the number of progenitors will be invariant, 
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Figure 9. The conditional mass function (CMF) for different descendant 
halo masses and at at different time-steps as indicated. The figure demon- 
strates that the transformation Au; = 1.3 X Au) properly compensates for 
the change in descendant mass Mq . 
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This equation can be used to find the transformation M — ► M nu- 
merically using the global fitting function of eq.fS] In fig.UOIwe 
show that this transformation can yield perfect matching between 
very different CMF's. The only limitation in the accuracy is the 
goodness of our global fit, which is used to find the mass transfor- 
mation above. As shown in fig. [10] the same mass transformation 
is suitable for a large range of time-steps, so the deviations of the 
CMF fit at small time-steps do not degrade the accuracy of the mass 
transformation. 

In general, one might use transformation of Au; and mass ac- 
cording to the specific case needed. Time transformation are much 
more easy to perform, and should be preferred for cases where the 
resulting time resolution is not problematic. Mass transformation 
can be done as a second step, to increase the accuracy of the statis- 
tics by small variations in mass. For exa mple, as tested above, using 
the results of the Millennium simulation dSpringel et alj2005h with 
a transformation Auj = 0.86Au; should give merger trees which 
are consistent with the cosmological model WMAP5 for a large 
range of halo masses. For this simple case, the time transformation 
is consistent with the naive EPS prediction and scales like erg. 



5.2 Generating Monte-Carlo trees 

Constructing Monte-Carlo merger trees is doable once an accu- 
rate knowledge of the CMF is given for all time-steps. There are 
various algorithms that were applied to the EPS CMF, and can 
be easily generalized for other CMF's l|Kauffmann & White|[l993l : 
iNeistein & Dekelll2008bl : IZhang et alj|2008l) . In this sense our fit- 
ting function for the CMF might be very useful for constructing 
Monte-Carlo trees that fit the results of iV-body simulations. How- 
ever, there are some limitations to this approach which we will em- 
phasize here below. 

As stated in the previous section, the quality of our global fit 
is poor for small time steps (Au; < 0.5). On the other hand, the 
time-step which is convenient for constructing Monte-Carlo trees 
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Figure 10. Scaling the mass in order to compensate for the change in 
dN/dM. We use eq. [l0]to transform the mass of progenitors from the 
wmap5 simulation (solid blue lines) into the scdml model (dashed red), 
both at the same time Auj. The mass transformation M; — > Mi derives 
from the global fit at a specific time-step, Auj = 1. The same transforma- 
tion is used to transform the actual progenitor masses for all the time steps 
used here. The dashed-dotted lines show the resulting CMF values, which 
are similar to the target case. The descendant halo mass is 10 12 Il~ 1 Mq. 



is much smaller, around Auj = 0.1. Larger time-steps yield many 
progenitors in each merger event, and a large uncertainty in the 
predicted merging time of the haloes. As a result, we do not have 
an accurate fitting function for the CMF, which is appropriate for 
generating Monte-Carlo trees. More than that, an accurate fit for 
the CMF at small Auj will not provide a solution, because the be- 
haviour of merger trees is highly non-markov at such small time- 
steps. This means that applying the CMF in consecutive time-steps 
without proper correlation between steps will generate large devia- 
tions in the merger histories (ND08a). 

A possible solution to these problems is to look for a new CMF 
at small time-steps that will reproduce the CMF from A-body sim- 
ulations at big time-steps. The CMF at small Auj does not have to 
match A-body simulations, it can only be tested by applying it on 
few consecutive time-steps. Such a methodology was introduced by 
ND08a and found to produce good results. However, the resulting 
trees are fully markov, differing from A-body simulations. As a re- 
sult of the above complications we think that this problem deserve 
more room than what is left here, and we postpone it to a future 
work. 



6 SUMMARY AND DISCUSSION 

We provide an improved description of the merger histories of 
dark-matter haloes as measured from A-body simulations. Using a 
suite of cosmological simulations, and a range of descendant halo 
masses at different redshifts, we provide a robust fitting function 
for the conditional mass function of progenitors, the CMF. This fit 
emphasizes the self-similar and universal natu re of the merger trees 
and it is more accurate than earlier attempts jLacev & Cole|| 19931 ; 
ISheth & Tormenll2002l : iMoreno et al]|2008l : ICole et al.ll2008l) . The 
improved fit is owing to the functional form used, which gener- 
alizes the EPS formula with a term that breaks the symmetry be- 
tween the natural time and mass variables. An improved accuracy 



is achieved by allowing the best-fit parameters to explicitly depend 
on the descendant halo mass. We note that it might be that other 
fitting functions will describe the data in a similar accuracy, and 
that our formula will deviate from the simulation results at different 
halo mass ranges. However, the range in halo mass probed in this 
work dominates the contribution to galaxy formation processes. 

The accurate fitting function for the CMF can serve as a ba- 
sis for comparing different simulations, or different algorithms for 
construction merger trees. I t can also be used for testi ng new mod- 
ifications to the theory (e.g. lMaggiore & Rio tto 2009) and for con- 
strain ing the time evolution of the halo mass function dTinker et al.l 
120081) . We discuss a specific application of our fitting function, that 
is to scale merger trees that are extracted from a given A-body sim- 
ulation into a different cosmological model, or different ranges of 
mass or redshift. We demonstrate that these scalings provide more 
accurate merger trees than other methods that utilize Monte-Carlo 
realizations. The scaling method conserves the correlation of the 
progenitor masses between time-steps, a non-Markovian effect that 
is hard to mimic in Monte-Carlo generated trees. It would be worth- 
while to study the scaling of substructures along similar lines. Sub- 
structures should be affected mostly by the mass ratio of the sub- 
halo to its host halo, and by the dynamical time-scale within the 
host. Once these quantities are conserved, substructures might be 
easily scaled as well. Such an analysis of substructures requires 
simulations of higher resolution than those used here. 

While most predictions of the EPS approximation are inac- 
curate, we find that the self-similarity in time proposed by EPS 
is valid to an accuracy of a few percents. We show explictly 
that merger histories are self-similar when using the time variable 
uj — S c /D(t), with 5 C — 1.68 and D(t) the cosmological lin- 
ear growth rate. This sel f-sim ilarity, which has b een noticed by 
iNeistein & Dekell j2008ah and [Genel et al.l J2008h in the Millen- 
nium simulation I Springel et alj2005b . is verified here for different 
cosmological models. We show that this self-similarity implies that 
uj can be used for relating merger trees from different cosmologies. 

In this work we used as default a standard halo definition 
based the FOF algorithm with a linking length corresponding to 
20% of the mean interparticle separation. We tested in compari- 
son merger trees for haloes defined using a spherical overdensity 
(SO) method (for m ore details on how these trees are defined see 
iMaccio et alj|2008h . and found that they are less self-similar in 
time, and are thus not as suitable for obtaining a universal fitting 
function. Indeed, th e FOF haloes are kn own to provide a better 
match to the theory dLacev&CoIell 19941) . and they are therefore 
more com monly used in the relevant literature (e.g. Sheth & Tor- 
men 200ffi lCole et alj2008t) . 

The construction of a merger tree by linking haloes from dif- 
ferent simulation snapshots can be done in different ways, and the 
results may depend on the method adopted. A comparison of our re- 
sults and other results from the litera ture indicates that this is not an 
issue of major concern. For example, ICole et al.l < f2008h does not re- 
port a noticable dif ference between their CMF results and those of 
iBower et al. which use different algorithms for construct- 

ing merger trees. On the othe r hand, our progenitor d efinition is 
very similar to the one used bv lSheth & Tormenl < |2002|) . so the im- 
provement in the fitting quality is not likely to be due to deviations 
in the simulation results. Finally, merger trees that were constructed 
by different methods may differ at small time steps, but these dif- 



3 The merger trees bui lt with the GIF simulatio n used bv lSheth & Tormer] 
120021) are described in lKauffmann et alj \ 19991) 
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ferences tend to diminish when large time steps a re considered, e.g . 
the Az > 0.5 steps used in our analysis (see e.g. iGenel et alJioOS. 
and their figure 8b). 

Our current analysis of a suite of cosmological simulations 
confirms the lognormal nature of the main-progenitor distribution 
as a function of the natur al mass variable a (M), as reported by 
iNeistein &"p ekel ( 20083) from the Millenium simulation. The lin- 
ear dependence of the parameters of the lognormal distribution on 
the natural time variable log Auj is much simpler than the time de- 
pendence of the CMF fit described in the first part of our paper. 
This motivates a theoretical search for the origin of the lognormal 
distribution of main progenitors. Both methods used here to quan- 
tify merger-trees as measured from iV-body simulations emphasize 
the limits of current theories in capturing the details of dark-matter 
evolution. 
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APPENDIX A: ADDITIONAL STATISTICS 
Al Self -similarity in time 

In this section we provide more details on the accuracy of the self- 
similarity in time, as discussed in section 13.11 In fig. I Al I we plot 
average main-progenitor histories for our set of simulations. The 
histories are plotted for descendant haloes identified at za = 0, 1, 2. 
Comparing histories for different zo we see small deviations with 
Zo, reaching ~ 20% at Acu = 5 (for the scdml & lcdm simu- 
lations). It should be kept in mind that cosmic variance is non- 
negligible in this plot. We also plot the EPS prediction f or the 
main-progenitor histories, as given by INeistein et al .1 1 |2006|) . The 
difference between the analytical EPS prediction and the results of 
TV-body simulations seem to change slightly wi th halo mass. Thes e 
results are sometimes different from the study o flZhaoetalJ ( l2008l) . 
We note that these authors compared median values from simula- 
tio ns against average val ues from the EPS formalism, as estimated 
by Ivan den Bosch! J2002h . This makes their comparison less accu- 
rate than what is done here. 

The full distribution of the main-progenitor, Pi (AS] So, Aui), 
is plotted in fig. lA2l for descendant haloes selected at different zo- 
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Figure Al. Average main-progenitor histories at different times. Smooth 
lines show main-progenitor histories for descendant halo Mo identified at 
Zo = 0, 1, 2 (zo = 2 was omitted for the high-mass bin due to a small 
statistical sample). Symbols and e rror-bars are the pred ictions of the EPS 
formalism following the formula of' lNeistein et alj |2006|) , and are indepen- 
dent on zo- The symbols shapes are diamonds, squares, triangles, circles 
and pluses for the scdml, scdm2, Icdm, wmap5, and Millennium simula- 
tions respectively. It is evident that self-similarity in time is valid for all the 
simulations. The deviations of the EPS formalism with respect to TV-body 
simulations are similar in all cases, although slight trends with cosmology 
and halo mass can be seen. 



Here as well, self-similarity in time is shown to be accurate with 
the exception of the scdml simulation at zo — 0. 

In figs, IA3l & IA4l we plot the mutual distribution of the two 
most massive progenitors. We show two dimensional histograms 
for descendant haloes selected at zo = 0, 1,2. The results for 
scdml & scdm2 simulations show accurate similarity for differ- 
ent zo • However, the Icdm and wmap5 show deviations a the level 
of 10-20% in mass. These deviations decrease at larger time steps, 
so they might be connected to non-markov effects at small time- 
steps, and their variation with redshift. In addition, our treatment 
of 'blacksplash' haloes may affect the results (see the discussion in 
sectionUTt. 



Figure A2. The full distribution of the main-progenitor mass at different 
times. In each panel we plot the results from one simulation as indicated, 
where the descendant halo mass M is between 10 12 and 10 13 /i _1 M , 
and it is selected at zo = 0, 1, 2, 3 (solid, dashed, dashed-dotted, dotted 
lines respectively). Main-progenitors are followed backwards in time until 
Auj = 0.5 & 2 (blue and red curves respectively). 
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Figure A3. The mutual distribution of the two most massive progenitors 
(Mi, M2). Each panel shows results from one simulation as indicated. The 
descendant halo Mo is identified at zo = 0, 1, 2 (plotted as solid, dashed, 
dashed-dotted lines respectively). The contour lines are plotted for 7 & 30% 
of the maximum histogram values. Descendant mass is between 10 12 and 
10 13 h _1 M Q , and Auj = 0.3 in all cases. 



A2 Scaling merger trees 

In this section we further examine the scaling of merger trees 
as it was applied in section IA2I We only show mutual distri- 
butions of the two most massive progenitors, (Mi,M2), as the 
main-pro genitor distribution is highly constrained by the CMF 
(see e.g. iNeistein et al. In fig. IA5I we examine the time- 

transformation which was used to generate fig. [8] The results of 
the mass scaling as applied in fig.[T0]are shown here in fig, IA6l 
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Figure A4. Same as fig.|A3]but with Aa; = 0.8. 



Figure A6. The mutual distribution of the two most massive progenitors 
(Mi, M2), for mass matching. Mass transformation for wmap5 was done 
in the same way as in fig.fjo] Au) = 0.3, M = 10 12 h~ 1 M Q . 
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Figure A5. The mutual distribution of the two most massive progenitors 
(Mi, M2), for Aw matching. The same Au) transformation is done as in 
fig. [8] The contour lines are plotted for 13 & 50% of the maximum his- 
togram values. Descendant mass is indicated in units of /i _1 Mq, Aui = 1. 



A3 A global fit to the CMF 

In fig. [AT] we provide more tests to our global fit, in a similar way 
to fig.© 



APPENDIX B: A GLOBAL FIT TO THE 
MAIN-PROGENITOR DISTRIBUTION 

In section[4]we claim that the lognormal distribution can be easily 
described for all time-steps and descendant masses using eqs.[5]& 
[6] This is shown in fig. IB II where full distributions of the main- 
progenitor mass are plotted. In tables IBTI & IB2I we summarize the 
parameters of the fits. 
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Figure A7. Similar to fig. [6] but for different time-steps as indicated. 



Universal Merger Histories 

Table Bl. The coefficients used for the global lognormal fit, eqs.[3]&[6] 
The standard deviation of In AS obeys the equation <j p = (ai log So + 
a 2 ) log Au> + C13 log So + 04, where a; are given below for each cosmol- 
ogy. For the Millennium simulation (MS) we used the fit given by ND08a, 
and transformed it to depend on So . 
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Table B2. The coefficients used for the global lognormal fit, eqs.[5]&[6] 
The mean of In AS obeys the equation fi p = (61 log So + 62) log + 
£>3 log So + 64, where 6; are given below for each cosmology. For the Mil- 
lennium simulation (MS) we used the fit given by ND()8a, and transformed 
it to depend on So ■ 
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Figure Bl. The full distribution of the main-progenitor for different cos- 
mologies and time-steps. Descendant haloes are selected at Zo = 1 with 
mass between 10 12 and 10 13 /i -1 Mq. Main-progenitor mass is followed 
backward in time until Alu = 0.5, 2, 4. Histograms show the simulation 
data, smooth lines are generated using our global fit. 



